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[57] ABSTRACT 

An apparatus for acquiring signals emitted by a fetus, 
identifying fetal heart beats and determining a fetal heart 
rate. Multiple sensor signals arc outputted by a passive fetal 
heart rate monitoring sensor. Multiple parallel nonlinear 
filters filter these multiple sensor signals to identify fetal 
heart beats in the signal data. A processor determines a fetal 
heart rate based on these identified fetal heart beats. The 
processor includes the use of a figure of merit weighting of 
heart rate estimates based on the identified heart beats from 
each filter for each signal. The fetal heart rate thus deter- 
mined is outputted to a display, storage, or communications 
channel. A method for enhanced fetal heart beat discrimi- 
nation includes acquiring signals from a fetus, identifying 
fetal heart beats from the signals by multiple parallel non- 
linear filtering, and determining a fetal heart rale based on 
the identified fetal heart beats. A figure of merit operation in 
this method provides for weighting a plurality of fetal heart 
rate estimates based on the identified fetal heart beats and 
selecting the highest ranking fetal heart rale estimate. 

32 Claims, 7 Drawing Sheets 
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PASSIVE FETAL HEART RATE 
MONITORING APPARATUS AND METHOD 
WITH ENHANCED FETAL HEART BEAT 
DISCRIMINATION 

ORIGIN OF TUB INVENTION 

The invention described herein was jointly made in the 
performance of work under a NASA contract and under a 
NASA Grant, and is subject to Section 305 of the National 
Aeronautics and Space Act of 1958, as amended, Public Law 
85-568 (72 Slat. 435; 42 USC 2457) and Public Law 96-517 
(35 U.S.C. 200 ct scq. ). In accordance with 35 USC 202, the 
grantee elected not to retain title. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The invention relates in general to the determination of 
the fetal heart rate of a fetus inside an expectant mother 
through the identification of fetal heart beats, and more 
specifically, to an improved method and apparatus for 
detecting fetal heart beats from the signals generated by a 
passive fetal heart rale monitoring sensor and determining a 
fetal heart rate based on the identified ictal heart beats. 

2. Description of the Related Art 

Various methodologies have been employed in the indus- 
try to determine fetal heart rales within a noise- contaminated 
signal from a fetal heart monitoring sensor. For example, a 
well known method in the industry is bandpass filtering 
which removes low frequency noise from heart beat signals. 
However, bandpass filtering fails to adjust for background 
noise frequencies which change with respect to time. A more 
sophisticated methodology which adapts to changing back- 
ground noise is described in U.S. Pal. No. 4,781,200 to 
Baker. The Baker patent discloses a portable, non-in vasivc 
system for monitoring fetal heart rate employing adaptive 
cancellation and spectral analysis. The Baker system is 
essentially an ambulatory system which uses adaptive inter- 
ference cancelling to isolate fetal heart beats and frequency 
domain analyses of multi-channel signals to extract useful 
information pertaining to fetal well-being, including fetal 
heart rate and absence of fetal movement. Another system, 
described in U.S. Pat. No. 5,209,237 to Rosenthal, also 
employs a similar noise cancellation technique, but uses 
multiple reference sensors to detect various background 
noises for adaptive cancellation and subsequent time-do- 
main, cross-spectral frequency analysis. While there arc 
certain benefits in using adaptive cancellation and spectral 
analysis, the use of interference cancellation to identify fetal 
heart beats is dependent on the quality of the reference 
signals. Noise not included in the reference signal may not 
be filtered from the signals emanating from the abdomen. 
While generally known sources of noise, such as the moth- 
er’s heart beat and breathing, are easily identified and 
filtered, extraneous noise from other organs or from the 
external environment may be missed. Therefore, the result- 
ing signal may not have been sufficiently filtered of noise 
and an inaccurate fetal heart rate may have been determined 
from erroneous fetal heart beat data. 

The use of a linear prediction methodology has also been 
used in the industry as the basis for identifying fetal heart 
beats during the determination of a fetal heart rate. Such an 
approach is described in Pretlow III, R. A. and Stoughton, J., 
“Signal Processing Methodologies for an Acoustic Fetal 
Heart Rate Monitor,” Master’s Thesis, April, 1991, Old 
Dominion University (hereafter “Pretlow ct al”). However, 
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such a methodology, while an improvement over adaptive 
cancellation, is still sensitive to noise not related to the fetal 
heart beat. Moreover, linear prediction has a significant 
deficiency inherent in the linear nature of its filter structure 
5 which prevents it from accurately modeling the nonlinearity 
of the background noise. The use of a two-layer, feed 
forward neural network has been suggested for use in 
identifying heart beats and determining heart rates. vSuch a 
method is discussed in detail in IIu, Y. II., Tompkins, W. J., 
to and Xuc, Q., “Artificial Neural Network for ECO Arrhyth- 
mia Monitoring,” Proceedings of IJCNN Inti. Joint Conf. on 
Neural Networks, 1992, vol. 2 (hereafter “IIu ct al”). 
However, the drawback of the methodology discussed by I lu 
ct al. is that it concentrates on eliminating noise from BCG 
15 data and docs not deal with supporting an acoustic-based 
monitoring apparatus. Indeed, the ECG environment differs 
significantly from the acoustic signal environment for fetal 
heart beats in terms of the type and amount of noise and the 
strength of the heart beat signal. The previously-mentioned 
?.o Pretlow ct al. article acknowledges the possibility of using 
a neural network to perform signature matching to identify 
a fetal heart beat, but concludes that neural network pro- 
cessing works best in offline waveform processing, such as 
for ECG analysis, and has not been found practical for real 
25 time fetal heart beat detection. 


SUMMARY OF’ THE INVENTION 

Accordingly, it is an object of the invention to overcome 
30 the previously described problems and deficiencies in exist- 
ing technology for identifying fetal heart beats and deter- 
mining an accurate fetal heart rate. 

In addition, it is an object of the invention to provide an 
35 improved apparatus and method for identifying fetal heart 
beats and determining fetal heart rates from signals gener- 
ated from an acoustic fetal heart monitoring sensor. 

It is a further object of the invention to provide a more 
accurate apparatus and method for identifying a fetal heart 
40 heat, through the simultaneous nonlinear filtering of a plu- 
rality of acoustic signals, generated by a passive fetal heart 
rate monitoring sensor, by the use of a Teager Energy 
Operator and/or a neural network nonlinear filter or filters, 
and determining a fetal heart rate based upon the identified 
45 fetal heart beats. 

In order to achieve the foregoing and other objects, in 
accordance with the purposes of the present invention as 
described therein, an apparatus for identifying fetal heart 
beats and determining a fetal heart rate comprises a passive 
50 fetal heart monitoring sensor producing a plurality of acous- 
tic signals, multiple nonlinear filters simultaneously filtering 
the plurality of sensor signals, and means for determining 
fetal heart rates. As embodied herein, there is a passive fetal 
heart rate monitoring sensor which acquires acoustic signals 
55 emitted from a fetus inside a body and outputs a plurality of 
sensor signals; a signal processing device, which receives, 
amplifies, filters, multiplexes, and digitizes the plurality of 
sensor signals, and which outputs a plurality of processed 
signals; multiple parallel nonlinear filters, which receive the 
60 plurality of processed signals, filter each of the plurality of 
processed signals to identify fetal heart beats, and output a 
plurality of functional values indicative of the presence of 
fetal heart beats for each of the plurality of processed 
signals; and fetal heart rate determining means, responsive 
65 to the plurality of functional values indicative of the pres- 
ence of fetal heart beats, for determining and outputting a 
fetal heart rate. 
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A method for identifying fetal heart beats and determining 
fetal heart rates comprises multiple parallel nonlinear filter- 
ing of a plurality of acoustic sensor signals generated from 
a passive fetal heart rate monitoring sensor and steps for 
determining a fetal heart rate based on the identified heart 
beats. As embodied herein, the method comprises a step for 
acquiring acoustic signals emitted from the fetus and out- 
putting a plurality of sensor signals, a step for receiving, 
amplifying, filtering, multiplexing, and digitizing the plu- 
rality of sensor signals to produce a plurality of processed 
signals, a step for multiple parallel nonlinear filtering of each 
of the plurality of processed signals to identify the presence 
of fetal heart beats and to output a plurality of functional 
values indicative of the presence of fetal heart beats for each 
of the plurality of processed signals, and a step for deter- 
mining and outputting the fetal heart rate based on the 
plurality of functional values indicative of the presence of 
fetal heart beats. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention and the objects achieved by it will 
be understood from the description of the preferred embodi- 
ments, taken in conjunction with the accompanying draw- 
ings of which: 

FIG. 1 is an overview of the hardware platform for an 
embodiment of an apparatus according to the present inven- 
tion. 

FIG. 2 is a block diagram of the pre-processing hardware, 
the signal processing system, and the computer for an 
embodiment of an apparatus according to the present inven- 
tion. 

FIG. 3 is a flow chart illustrating the multiple parallel 
nonlinear filtering for the multiple sensor signals and the 
determination of a fetal heart rate for an embodiment of a 
method according to the present invention. 

FIG. 4 is a graph depicting the time domain acoustic fetal 
heart signal and the Teager Energy Operator response to the 
acoustic signal according to an embodiment of the present 
invention. 

FIG. 5 is a schematic view illustrating a neural network 
according to an embodiment of the present invention. 

FIG. 6 is a schematic view illustrating a particular unit in 
the neural network according to an embodiment the present 
invention. 

FIG. 7 is a graph showing the time domain acoustic fetal 
heart signal and the neural network response to the signal 
according to an embodiment of the present invention. 

FIG. 8 is a graph depicting the fetal heart rate tracing 
obtained from acoustic monitoring of a patient with a 
gestation age of 39 weeks. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

Reference will now be made in detail to the present 
preferred embodiments of the present invention, examples 
of which are illustrated in the accompanying drawings, 
wherein like reference numerals refer to like elements 
throughout. 

FIG. 1 is an overview of the hardware platform for an 
apparatus according to an embodiment of the present inven- 
tion. In FIG. 1, an ambulatory, non-invasive, passive, fetal 
heart rate monitoring sensor 1 is placed on the abdomen of 
an expectant mother. Multiple sensor elements inside sensor 
1 receive the signals emitted by a fetus inside the expectant 
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mother and output a plurality of corresponding sensor sig- 
nals 2, commensurate with the number of sensor elements 
inside sensor 1. Such ambulatory, non-invasive, sensor ele- 
ments are well known in the industry and an example of this 
5 type of sensor 1 is described in U.S. Pat. No. 5,140,992 to 
Zucherwar et al. entitled “PASSIVE FETAL MONITOR- 
ING SENSOR.” The plurality of sensor signals 2 are 
received by a preprocessing electronic unit 3 which ampli- 
fies, filters, and multiplexes the signals in order to prepare 
10 and organize the sensor signals 2 for further processing. The 
preprocessed signals outputted by the preprocessing elec- 
tronic unit 3 are received by the signal processing system 4 
which demultiplexes the preprocessed signals back into 
multiple sensor signals and digitizes the preprocessed mul- 
15 tiple sensor signals for further analysis by an improved 
method of identifying fetal heart beats and determining a 
fetal heart rate. In particular, multiple parallel nonlinear 
filters are used in achieving significantly more accurate 
results for the identification of fetal heart beats. Moreover, 
20 the determination of multiple fetal heart rate estimates, one 
estimate from the result of each nonlinear filtering of each 
sensor signal, and the use of a figure of merit method to 
select the best fetal heart rate estimate ensures the most 
accurate result from this system. The signal processing 
25 system 4 is programmed to carry out these functions of heart 
beat identification and heart rate determination. 

As to the signal processing system 4, one embodiment of 
the present invention has the signal processing system 4 
inside a personal computer (PC) 5. The PC 5 is well known 
30 in the industry and, as shown in FIG. 2, is comprised of a 
central processing unit (CPU) 12, memory 13, and fixed disk 
storage 14; and further, has attached a keyboard 15 and a 
display 16. The PC 5 is used to carry out input from and 
output to various devices, such as the keyboard 15, the 
35 display 16, or a communications channel (not shown). The 
signal processing system 4 is not limited to being physically 
located inside a PC. For example, the preprocessing elec- 
tronic unit 3 and the signal processing system 4 can be 
combined in a signal processing device separate from the 
40 PC. In addition, the embodiments of the present invention 
are not limited to the identification of fetal heart beats and 
determination of fetal heart rates by a programmed signal 
processing system 4. For example, a general purpose digital 
computer, such as a generic version of PC 5, may be 
45 programmed to provide the same features. After a fetal heart 
rate is determined, the result may be displayed on display 16 
or stored in fixed storage 14 in the PC, for example, when 
appropriate commands are entered via the keyboard 15. 

FIG. 2 is a block diagram detailing the preprocessing 
50 electronic unit 3, the signal processing system 4, and the PC 
5. The preprocessing electronic unit 3 comprises amplifiers 
6 which receive the plurality of sensor signals 2. The sensor 
signals 2 outputted by the passive fetal heart rate monitoring 
sensor 1 is weak due to the inherent weakness of a fetal heart 
55 beat. One purpose of the amplifiers 6 is to boost the initial 
sensor signals 2. Another purpose of the amplifiers 6 is to 
amplify the sensor signals 2 to match the range of subse- 
quent components, e.g. an analog-to-digital converter. The 
amplifiers 6 respectively are connected to corresponding 
60 filters 7. The sensor signals 2 are bandpass filtered by 
bandpass filters 7 in order to eliminate competing signals 
that are larger than the fetal heart beat signal. By focusing 
the bandwidth on the relevant range of a fetal heart beat, the 
analog-to-digital converter and other subsequent processing 
65 will not be saturated with non-essential signals. The band- 
pass filters 7 are connected to an analog time-division 
multiplexor (M/X) 8. The multiplexor rotates among the 
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filtered signals and time-division multiplexes the filtered 
signals to allow simultaneous sampling of all sensor signals. 
The time slice lor each signal is very small (c.g., 0.5 ms) to 
accommodate the monitoring of dynamic, real-time changes 
in fetal heart rates. The multiplexor 8 is connected to an 
analog-to-digital converter 9 in the signal processing system 
4. In one embodiment, the signal processing system is a 
digital signal processing (DSP) card housed inside PC 5. The 
signal processing system 4 has a central processing unit 
(CPU) 10 and memory 11 in addition to signal processing 
electronics (c.g., analog-to-digital converter 9). However, 
the embodiments of the present invention arc not limited to 
having a signal processing system 4 as a DSP card residing 
inside a PC. For example, the preprocessing electronic unit 

3 may be combined with the features and functions of the 
signal processing system 4 in a separate signal processing 
device outside the PC. As a further example, a general 
purpose digital computer, such as a generic version of PC 5, 
may be programmed to provide the same features as the 
signal processing system 4. The prcproccsscd signal output- 
ted from the preprocessing electronics unit 3 is received by 
the signal processing system 4. The signal processing system 

4 first converts the analog prcproccsscd signal into digital 
words, at a specified sampling rate, c.g., 500 11/ or wave- 
form samples per second. Processing performed by the DSP 
CPU 10 receives the digital words via data stored in the DSP 
memory 11 and demultiplexes the digitized signal so that 
further processing may address time slices of digitized 
acoustic data for each of the multiple sensor signals 17 (as 
shown in FIG. 3). 

FIG. 3 is a flow chart showing further processing by the 
DSP CPU 10 of the signal processing system 4. The digi- 
tized signal is further filtered (c.g., at a 500 samples per 
second rate) with a linear phase FIR digital filter 18 (“FIR 
LPF”). The use of digital fillers operating on the digitized 
signal reduces the requirements on the analog bandpass 
filters 7 used in the hardware of the preprocessing electron- 
ics unit 3 and allows for greater flexibility since these filters 
arc very easily changed or tailored to any reasonable desired 
response. The linear phase property of the FIR digital filter 
also insures minimal waveform distortion. The filtered digi- 
tal signal may then be further decimated to a lower sampling 
rate (c.g., 250 samples per second). Such decimation to a 
lower sampling rate may be desirable since it reduces the 
computations required in subsequent steps of processing. 

After the digital filtering of the multiple sensor signal data 
17 by filters 18, the signals arc further analyzed to identify 
fetal heart beats and a resultant fetal heart rate is determined 
from the identified fetal heart beats. Multiple parallel non- 
linear filters arc used to identify fetal heart beats from the 
digitized acoustic sensor data. FIG. 3 shows two such 
nonlinear filters, i.c. a Teager Energy Operator 19 and a 
neural network 20. These nonlinear fillers operate in parallel 
and provide redundant identification of fetal heart beats for 
each of the plurality of sensor signals 17. Such parallel 
redundant processing improves the accuracy of identifying 
the fetal heart beats. While these benefits arc created by the 
use of multiple parallel nonlinear filters, the embodiments of 
the invention arc not limited to multiple parallel nonlinear 
filters. For example, singular nonlinear filtering, by the 
Teager Energy Operator 19 or by the neural network 20, of 
the plurality of sensor signals 17 is also possible to yield 
improved identification of fetal heart beats. 

The particular embodiment in FIG. 3 shows a first non- 
linear filtering by a Teager Energy Operator 19. The Teager 
Energy Operator 19 highlights regions of high energy (i.c., 
the heart beats). This nonlinear filler is known in the industry 
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and is described in Kaiser, James F., “On a Simple Algo- 
rithm to Calculate the ‘Energy’ of a Signal,” ICASSP-90, 
1990 , vol. 1 , pp. 381 -84 (hereafter “Kaiser”). The Teager 
Energy Operator is computed using 
5 

;y(/i) -x\n- 1) jc(m)*jc<« 2), 

where x(n) is the acoustic signal and y(n) is the Teager 
energy signal. This nonlinear operation results in a large 
10 output in regions which contain high amplitude and rapidly 
varying signals. The Teager signal is further digitally 
smoothed using a “box-car” average filter, well known in the 
industry. FIG. 4 depicts an example of an acoustic signal 28 
and the Teager Energy Operator filtered response signal 29. 
15 The acoustic signal 28 shown in FIG. 4 is the amplitude of 
the acoustic signal as time passes. The Teager Energy 
Operator filtered response signal 29 is the amplitude of the 
Teager Energy signal as time passes. 

The other nonlinear filter shown in FIG. 3 is a neural 
20 network. In one embodiment of the present invention, the 
neural network 20 is a three-layer feed- forward neural 
network as shown schematically in FIG, 5. The three layers 
of the network include an input layer 30, a hidden layer 31, 
and an output layer 32. For this nonlinear filter, the data is 
25 first normalized by a root mean square (RMS) value to 
insure that the subsequent processing receives data of uni- 
form amplitude. In other words, the RMS value of each 
sampled signal may be used to control or adjust the input 
range of the analog-to- digital convertor to match the 
30 measured signal level. Such a normalization method helps to 
maintain the best possible signal-to-noisc ratio for the signal 
data. The data is then processed by the neural network 20 
which nonlincarly “filters” the data using twenty digital 
words 33 (corresponding to the current input point plus 
35 nineteen previous input points 36) to compute an output 
point 36 for the current time frame, using an equation: 


yin) f{x{n), *(« 1 ), . . . *(« 19 )}, 

40 for each index in the frame of data, 
where 

y(n) is the current neural network output data 36 
f represents the neural network 
45 x(n) is the current input data point 36 

In addition, the neural network 20 comprises thirty hidden 
layer units 34 and one output layer unit 35. FIG. 6 is a 
schematic depiction of an arbitrary unit (“Unit X”) corre- 
sponding to a hidden layer unit 34 or to the output layer unit 
50 35. With unit X corresponding to a particular hidden layer 
unit 34, the plurality of inputs 37 (X r -X„) correspond to the 
plurality of digital words 33 from the input layer 30. Each 
input 37 is weighted by a weight 38 (W r W n ), summed with 
a bias term 39, and processed by a sigmoidal nonlinearity 40. 
55 An output value 41 is outputted to the output layer unit 35. 
With unit X corresponding to the output layer unit 35, each 
of the plurality of inputs 37 (X r -X„) corresponds to the 
outputted value from each of the hidden layer units 34. An 
output value 41, corresponding to the neural network output 
60 value 36, is outputted by the output layer unit 35. Thus, the 
neural network 20 is fully interconnected with a weight from 
each input layer unit 30 to each hidden layer unit 34, and 
from each hidden layer unit 34 to the output layer unit 35, 
and a bias term 39 for each unit. Bipolar sigmoidal nonl in- 
65 caritics 40 arc used for all hidden layer units 34 and the 
output layer unit 35. In equation form, the neural network 
can be described as follows: 
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The input to each hidden layer unit 34 is of the form: 
19 

8i= ‘WijxW + bi 
7=0 

for 0^29 
where 

g^input to \ th hidden unit 
W 0 -=weight from input j to hidden unit i 
x(j)=j"‘ input 
b—bias for hidden unit i 

The output of each hidden layer unit 34 is of the form: 


where 

g=input to \ th hidden unit 
h—output of \ th hidden unit 

The input to the final output layer unit 35 is of the form: 

29 

0= S whjhj + b 
7=0 

where 

0=input to output unit 

wh,— weight from f l hidden unit to output unit 35 
h^ output of jth hidden unit 
b=bias for output unit 35 

The final output, i.e., from the output layer unit 35, is 
given by: 

P 1 + e~° 

where 

p=output of output unit 35 
0=input of output unit 
The neural network 20 was trained with a back propaga- 40 
tion method using hand labeled training data, wherein 
regions of fetal heart beats and between-beat intervals were 
previously identified. Such data were obtained from clinical 
tests. Backpropagation is a well known method in the 
industry used to train a neural network, as described in 45 
Lippmann, R. P., “An Introduction to Computing with 
Neural Nets,” IEEE ASSP magazine, April 1987, pp. 4-22 
(hereafter “Lippman”). The values for the weights and bias 
terms were determined from this training procedure. After 
the completion of the training, the neural network 20 oper- 50 
ates using the learned weights and biases. In essence, the 
neural network 20 acts as a nonlinear filter trained to indicate 
the presence of a fetal heart beat. FIG. 7 is a graph showing 
an example of a fetal heart tone as represented by an acoustic 
signal 42 and the resulting, nonlinearly filtered time-domain 55 
output waveform 43 from the neural network 20. The 
acoustic signal 42 shown in FIG. 7 is the amplitude of the 
acoustic signal as time passes. The neural network output 
waveform 43 is the amplitude of the neural network output 
as time passes. 60 

As regards the neural network 20, although a layered, 
interconnected, feed forward neural network is caused to 
learn by back propagation, the embodiments of the present 
invention are not limited to this choice of a neural network 
nor to this choice of training. For example, an interlinked 65 
feed forward neural network also having feed back connec- 
tions can be used. Such a feed forward and feed back neural 
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network can be trained by a method proposed by J. J. 
Hopfield or another by G. E. Hinton et al. These learning 
methods are described in Hopfield, J. J., “Neural network 
and physical systems with emergent collective computa- 
5 tional abilities,” Proc. Natl. Acad. Sci. U.S.A., 1982, vol. 79; 
Hopfield, J. K., “Neurons with graded response have col- 
lective computational properties like those of two-state 
neurons,” Proc. Natl. Acad. Sci. U.S.A., 1984, vol. 81; 
Ackley, D. H., Hinton, G. E„ and Senjnowski, T. J., “A 
10 Learning Algorithm for Boltzmann Machines,” Cognitive 
Sci., 1985, 9. Other variations of the neural network are also 
possible. For example, neural networks having a different 
number of hidden units, output units, or layers may also be 
15 trained to identify fetal heart beats. 

The resulting signal from each of the multiple parallel 
nonlinear filtering is then analyzed with autocorrelation 
processing 21 (as shown in FIG. 3) as a first step in 
determining a fetal heart rate estimate 23. The determination 
20 of a fetal heart rate estimate for each of the plurality of 
results from the multiple parallel nonlinear filtering of each 
of the plurality of sensor signals 17 will increase the chances 
of discovering the most accurate fetal heart rate estimate. 

Peak evaluation processing 22, as the next step as shown 
25 in FIG. 3, identifies peaks, selects the highest peaks, esti- 
mates a fetal heart rate based on the selected peaks in the 
autocorrelation, and determines a figure of merit 24 for the 
fetal heart rate estimate 23. A peak is defined as a point in 
the autocorrelation which is larger than immediately adja- 
30 cent points. A peak is selected for fetal heart rate estimation 
if the peak is larger than a threshold number of adjacent 
points. A heart rate estimate is found by identifying and 
using the selected peaks in the autocorrelation over a search 
35 range appropriate for fetal heart rates (typically 80 to 200 
beats per minute). The heart rate estimates 23 are determined 
in this manner for each of the multiple sensor channels. The 
heart rate estimates 23 are then ranked using a figure of merit 
24, which is derived from continuity constraints and a 
measure of periodicity in the waveform, as described below. 
The highest ranking (i.e., having the highest figure of merit) 
estimate is then selected in the next step by the merit 
comparison processing 25 as the heart rate 26 for the current 
frame of data. The figure of merit 27 is the figure of merit 
value corresponding to the selected fetal heart rate 26. This 
procedure is repeated for each frame of data. The time 
between frames for the entire processing is very small (e.g., 
0.5 sec.) in order to properly monitor dynamic real-time 
changes in fetal heart rates. 

The figure of merit 24 is determined as follows. First, one 
buffer is used to store a small number (e.g., 5) of the 
previously selected values of fetal heart rate 26. Another 
buffer is used to store a corresponding small number (e.g., 
5) of the previously corresponding values of figure of merit. 
A small number of previously calculated fetal heart rates and 
corresponding merit values is used to ensure the best 
dynamic response to fetal heart rate changes. A calculation 
is made using the following variables and equations: 

For i ih rate estimate, 

Merit = Peak(i) * ^ 1.0 - Ra t eDe^ ) * (Rate ® “ Av g- Rate > 
where 

Merit=Figure of merit for the i th fetal heart rate estimate. 
Peak(i)=The normalized autocorrelation for the \ th fetal 
heart rate estimate. 
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1 m 

Avg. Merit*- — —7 — X Mcrit(i) 
wh- 1 j'O 

where 5 

m is a small number (c.g., 4, which will provide for five 
samples when i ranges from 0 to 4), and 
Mcrit(i)=prcvious merit computations 27 
Rate Dcv.=Thc factor which controls the contribution of 
previous estimates (typically— 30). to 

Rate(i)=fctal heart rate estimate 23 

m 

X Merit (0 ■ Ratc(f') 

Avg. Rate 

6 m k 

X Mcrit(0 15 

i 0 

where 

m is a small number (c.g., 4, which will provide for five 
samples when i ranges from 0 to 4); 20 

Mcrit(i)=prcvious merit computations 27; and 
Ratc(i)=prcvious selected rates 26. 

This Avg. Rate is the “weighted” average of a small 
number m of previously selected fetal heart rates 26 having 
the highest ranking merits 27 when these rates 26 were ' " 
selected by the merit comparison operation 25. 

The heart rale estimate 23 with the highest figure of merit 
24 is selected as the current and most accurate fetal heart rale 
26 during merit comparison 25. If the figure of merit 27 for 
this calculation is above an empirically determined threshold 3 
(typically about 35% of the maximum value), the resultant 
calculation is assumed to be valid. When and if the output 
heart rate 26 is determined to be valid and is to be displayed, 
it is displayed as a “correct” reading. If the merit 27 is lower 
than the empirical threshold value (i.c., an invalid rate), the 3 
rate 26 is still displayed, but is color coded in the computer 
display to indicate a possible error or “dropout.” 

Alternatively, the outputted fetal heart rates 26 may be 
saved in fixed disk storage on the PC or transmitted to a 
communications channel, such as a modem. An example of 4 
an output tracing of fetal heart rates from a device for an 
embodiment according to the present invention, obtained 
from a patient with a gestation age of 39 weeks, is shown in 
FIG. 8. The y-axis is measured in beats per minute versus 
time as measured along the X-axis. Hach gradation along the 4 
X-axis equals approximately 20 seconds, or three gradations 
for about a minute. 

Numerous modifications and adaptations of the present 
invention will be apparent to those skilled in the art. Thus, 
the following claims and their equivalents arc intended to 50 
cover all such modifications and adaptations which fall 
within the true spirit and scope of the present invention. 

What is claimed is: 

1. A passive fetal heart rate monitoring apparatus com- 
prising: 33 

a passive fetal heart rate monitoring sensor having a 
plurality of sensor elements lor acquiring acoustic 
signals emitted from a fetus inside a body and output- 
ting a plurality of sensor signals in response to the 6() 
acoustic signals; 

a signal processing device, including receivers, amplifi- 
ers, and filters for outputting a plurality of processed 
signals in response to the sensor signals; 
means for detecting fetal heart beats comprising a neural 65 
network, responsive to the plurality of processed sig- 
nals, for identifying regions of fetal heart beats on each 
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signal, and for outputting a functional value indicative 
of the presence of fetal heart beats for each processed 
signal; 

fetal heart rate determining means, responsive to the 
plurality of functional values indicative of the presence 
of fetal heart beats, for determining and outputting a 
fetal heart rate. 

2. An apparatus as in claim 1, wherein said neural network 
comprises an input layer, a hidden layer, and an output layer, 
each layer comprising one or more nodes, a first node in the 
input layer being connected to the current input processed 
signal, and all other nodes in the input layer being connected 
to a different previous input point corresponding to a pre- 
vious input processed signal, each node in the input layer 
being connected to each node in the hidden layer and each 
node in the hidden layer being connected to each node in the 
output layer, the output layer outputting a functional value 
indicative of fetal heart beats, each connection between 
nodes having an associated weight; and 

a training means for determining the weight for each said 
connection between nodes of the neural network. 

3. An apparatus as in claim 2 wherein the training means 
comprises means for adjusting each weight by use of back 
propagation. 

4. An apparatus as in claim 1, wherein said means for 
detecting fetal heart beats further comprises a plurality of 
neural networks connected in parallel to said neural network 
for providing redundant identification of the fetal heart beats 
from the signal, and outputting a plurality of functional 
values indicative of the presence of fetal heart beats for each 
processed signal. 

5. An apparatus as in claim 1, wherein said means for 
detecting fetal heart beats further comprises a Teager Energy 
Operator connected in parallel to said neural network for 
providing redundant identification of the fetal heart beats for 
each received processed signal by outputting a plurality of 
functional values indicative of the presence of fetal heart 
beats for each processed signal, said Teager Energy Operator 
is defined by the equation: 


y(w) x~2*(n 1) *(«)*(«• -2) 

where x(n) is the processed signal and y(n) is the functional 
value. 

6. An apparatus as in claim 1 wherein the fetal heart rate 
determining means comprises: 

means for determining and outputting a fetal heart rate 
estimate for each of the received functional values; 

means for ranking the estimated fetal heart rates by 
determining a figure of merit value for each received 
fetal heart rate estimate; and 

merit comparison processing means for comparing the 
received figure of merit values for each received esti- 
mated fetal heart rale and outputting the estimated fetal 
heart rate which has the highest figure of merit. 

7, An apparatus as in claim 6, wherein said means for 
determining and outputting a fetal heart rate estimate com- 
prises: 

an autocorrclator for forming an autocorrelation function 
responsive to each of the received plurality of func- 
tional values, and 

means for delecting peaks for each received autocorrela- 
tion function over a search range appropriate for fetal 
heart rates, converting the peaks into an estimated fetal 
heart rate for each autocorrelation function, and out- 
putting each of the plurality of estimated heart rates. 



5,524,631 


11 

8. An apparatus as in claim 1 further comprising display 
means for receiving and displaying the estimated fetal heart 
rate with the highest figure of merit. 

9. An apparatus as in claim 1 wherein said signal pro- 
cessing device further comprises an analog to digital con- 5 
verter for digitizing the plurality of sensor signals. 

10. A passive fetal heart rate monitoring apparatus com- 
prising; 

a passive fetal heart rate monitoring sensor having a 
plurality of sensor elements for acquiring acoustic 10 
signals emitted from a fetus inside a body and output- 
ting a plurality of sensor signals in response to the 
acoustic signals; 

a signal processing device, including receivers, amplifi- 
ers, and filters for outputting a plurality of processed 15 
signals in response to the sensor signals; 

means for detecting fetal heart beats comprising a Teager 
Energy Operator, responsive to the plurality of pro- 
cessed signals, for identifying regions of fetal heart 
beats on each signal, and for outputting a functional 20 
value indicative of the presence of fetal heart beats for 
each processed signal; 

fetal heart rate determining means, responsive to the 
plurality of functional values indicative of the presence 
of fetal heart beats, for determining and outputting a 25 
fetal heart rate. 

11. An apparatus as in claim 10 wherein said Teager 
Energy Operator is defined by the equation: 

30 

y(n)~xA2* («-l )-*(n)*(n-2) 

where x(n) is the processed signal and y(n) is the functional 
value. 

12. An apparatus as in claim 10 wherein said means for 35 
detecting fetal heart beats further comprises Teager Energy 
Operators connected in parallel to said Teager Energy 
Operator for providing redundant identification of the fetal 
heart beats for each received processed signal by outputting 

a plurality of functional values indicative of the presence of 40 
fetal heart beats for each processed signal. 

13. An apparatus as in claim 10 wherein the fetal heart rate 
determining means comprises: 

means for determining and outputting a fetal heart rate 
estimate for each of the received functional values; 45 

means for ranking the estimated fetal heart rates by 
determining a figure of merit value for each received 
fetal heart rate estimate; and 

merit comparison processing means for comparing the 5Q 
received figure of merit values for each received esti- 
mated fetal heart rate and outputting the estimated fetal 
heart rate which has the highest figure of merit. 

14. An apparatus as in claim 13, wherein said means for 
determining and outputting a fetal heart rate estimate com- 55 
prises: 

an autocorrelator for forming an autocorrelation function 
responsive to each of the received plurality of func- 
tional values, and 

means for detecting peaks for each received autocorrela- 60 
tion function over a search range appropriate for fetal 
heart rates, converting the peaks into an estimated fetal 
heart rate for each autocorrelation function, and out- 
putting each of the plurality of estimated heart rates. 

15. An apparatus as in claim 10 further comprising display 65 
means for receiving and displaying the estimated fetal heart 
rate with the highest figure of merit. 
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16. An apparatus as in claim 10 wherein said signal 
processing device further comprises an analog to digital 
converter for digitizing the plurality of sensor signals. 

17. A method for passive fetal heart rate monitoring of a 
fetus in a body comprising: 

acquiring acoustic signals emitted from a fetus inside a 
body and outputting a plurality of sensor signals in 
response to the acoustic signals; 

receiving, amplifying, and filtering outputting a plurality 
of processed signals in response to the sensor signals; 

detecting fetal heart beats using a neural network, respon- 
sive to the plurality of processed signals, identifying 
regions of fetal heart beats on each signal, and output- 
ting a functional value indicative of the presence of 
fetal heart beats for each processed signal; 

determining and outputting the fetal heart rate in response 
to the plurality of functional values indicative of the 
presence of fetal heart beats. 

18. The method of claim 17, wherein said neural network 
comprises an input layer, a hidden layer, and an output layer, 
each layer comprising one or more nodes, a first node in the 
input layer being connected to the current input processed 
signal, and all other nodes in the input layer being connected 
to a different previous input point corresponding to a pre- 
vious input processed signal, each node in the input layer 
being connected to each node in the hidden layer and each 
node in the hidden layer being connected to each node in the 
output layer, the output layer outputting a functional value 
indicative of fetal heart beats, each connection between 
nodes having an associated weight; and 

a training means for determining the weight for each said 
connection between nodes of the neural network. 

19. The method of claim 18 wherein the training means 
comprises means for adjusting each weight by use of back 
propagation. 

20. The method of claim 17, wherein the step of detecting 
fetal heart beats further comprises the steps of providing 
redundant identification of the fetal heart beats from each 
signal using a plurality of neural networks connected in 
parallel to said neural network and responsive to each signal, 
and outputting a plurality of functional values indicative of 
the presence of fetal heart beats for each processed signal. 

21. The method of claim 17, wherein the step of detecting 
fetal heart beats further comprises the steps of providing 
redundant identification of the fetal heart beats for each 
received processed signal using a Teager Energy Operator 
connected in parallel to said neural network and responsive 
to each signal, and outputting a plurality of functional values 
indicative of the presence of fetal heart beats for each 
processed signal, said Teager Energy Operator is defined by 
the equation: 

y(n)=xA2*(n-l)-x(n)*(n-2) 

where x(n) is the processed signal and y(n) is the functional 
value. 

22. The method of claim 17 wherein the fetal heart rate 
determining means comprises: 

determining and outputting a fetal heart rate estimate for 
each of the received functional values; 

ranking the estimated fetal heart rates by determining a 
figure of merit value for each received fetal heart rate 
estimate; and 

comparing the received figure of merit values for each 
received estimated fetal heart rate and outputting the 
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estimated fetal heart rate which has the highest figure of 
merit. 

23. The method of claim 22, wherein the step of deter- 
mining and outputting a fetal heart rate estimate comprises: 

forming an autocorrelation function responsive to each of 5 
the received plurality of functional values, 

detecting peaks for each received autocorrelation function 
over a search range appropriate for fetal heart rates, 
converting the peaks into an estimated fetal heart rate 
for each autocorrelation function, and outputting each 10 
of the plurality of estimated heart rates. 

24. The method of claim 17 further comprising the steps 
of receiving and displaying the estimated fetal heart rate 
with the highest figure of merit. 

25. The method of claim 17 wherein the step of processing 1 5 
the signals further comprises the step of digitizing the 
plurality of sensor signals. 

26. A passive fetal heart rate monitoring method com- 
prising the steps of: 

acquiring acoustic signals emitted from a fetus inside a 
body and outputting a plurality of sensor signals in 
response to the acoustic signals; 

processing the signals using a signal processing device 
having receivers, amplifiers, and filters, and outputting 25 
a plurality of processed signals in response to the 
sensor signals; 

detecting fetal heart beats using a Teager Energy Operator 
responsive to the plurality of processed signals, iden- 
tifying regions of fetal heart beats on each signal, and 30 
outputting a functional value indicative of the presence 
of fetal heart beats for each processed signal; 

determining and outputting a fetal heart rate in response 
to the plurality of functional values indicative of the 
presence of fetal heart beats. 35 

27. The method of claim 26 wherein said Teager Energy 
Operator is defined by the equation: 
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where x(n)is the processed signal and y(n)is the functional 
value. 

28. The method of claim 26 wherein the step of delecting 
fetal heart beats further comprises the steps of providing 
redundant identification of the fetal heart beats for each 
received processed signal by connecting a plurality of Tea- 
ger Energy Operators in parallel to said Teager Energy 
Operator and outputting a plurality of functional values 
indicative of the presence of fetal heart beats for each 
processed signal. 

29. The method of claim 26 wherein the fetal heart rate 
determining means comprises: 

determining and outputting a fetal heart rate estimate for 
each of the received functional values; 

ranking the estimated fetal heart rates by determining a 
figure of merit value for each received fetal heart rate 
estimate; and 

comparing the received figure of merit values for each 
received estimated fetal heart rate and outputting the 
estimated fetal heart rate which has the highest figure of 
merit. 

30. The method of claim 29, wherein the step of deter- 
mining and outputting a fetal heart rate estimate comprises: 

forming an autocorrelation function responsive to each of 
the received plurality of functional values, 

detecting peaks for each received autocorrelation function 
over a search range appropriate for fetal heart rales, 
converting the peaks into an estimated fetal heart rate 
for each autocorrelation function, and outputting each 
of the plurality of estimated heart rates. 

31. The method of claim 26 further comprising the steps 
of receiving and displaying the estimated fetal heart rale 
with the highest figure of merit. 

32. The method of claim 26 wherein the step of processing 
the signals further comprises the step of digitizing the 
plurality of sensor signals. 


y{n) -xA2*(n - 1 ) -Jt(n)*(n- - 2 ) 





